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Time-stationary solutions to the Vlasov-Poisson equation for ion holes 
and double layers are examined along with particle simulations which pertain to 
recent observations of small amplitude <?$/ 7^—1 electric field structures on 
auroral field lines Both the time-stationary analysis and the simulations sug- 
gest that double layers evolve from holes in ion phase space when their ampli- 
tude reaches e<f>/ T e ~1. Multiple small amplitude double layers which have 
been seen in long simulation systems and are seen to propagate past the space- 
craft may account for the acceleration of plasma sheet electrons to produce the 
discrete aurora. 


We have examined time-stationary solutions to the Vlasov-Poisson equation for ion holes 
and double layers which pertain to particle simulations we have performed and recent observa- 
tions in space [Temertn et al , 1982] In the simulations, small amplitude e<f>/ T e ~1 double 
layers appear to evolve from holes in ion phase space when a current carried by electrons drift- 
ing at or near their thermal speed is applied to a long system — 1000 Debye lengths In longer 
current-driven simulation systems, multiple e<f>/ T e ~1 double layers evolve [Sato and Okuda, 
1981], In contrast, when a potential is applied across a simulation system, whichever e<f>/ T e is 
applied appears as an isolated double layer which forms at the boundary and propagates into the 
system [Hubbard and Joyce, 1979; Kmdel et al 1981] 

Figure 1 shows an example of small amplitude e$j T e ~ 1 double layers and solitary 
waves observed on the S3-3 satellite at 5000 km over the auroral zone There are two types of 
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aperiodic structures in the parallel electric field component which appear to be independent of 
the ion cyclotron wave period in the perpendicular electric field component symmetric solitary 
wave pulses with no net potential jump and asymmetric pulses which we identify as small 
amplitude e<j>/ T e double layers These structures are observed on the millisecond times- 
cale to propagate past the spacecraft at a velocity greater than 50 km/s For a typical example, 
the inferred velocity was 80 km/s and the pulsewidth was 4 ms This yields a scale length of 
320 m, which is 43 Debye lengths assuming a 10 eV, 10 cm -3 thermal electron population. 
This scale length for the double layer thickness is consistent with the current-driven simula- 
tions, as is the spacing of one to a few thousand Debye lengths 

The propagation speed, a lower limit of 50-100 km/s, is comparable to both the energetic 
ion beam speed and the ion acoustic speed in the region where the double layers are observed. 
These are regions of upward-directed ion beams and downward-accelerated electrons measured 
by the particle detectors on S3-3. Temerm et a I [1982] suggest that it is possible to account for 
a substantial portion of the potential drop necessary to explain the observed particle distribu- 
tions in the cumulative potential drop of the multiple double layers along an auroral field line. 
It is not possible to determine the direction of propagation from the double probe measure- 
ment, as explained in Figure 3 of Temerm et al. [1982]. However, we expect the solitary 
wave/double layer structures to be tied to the ion beam [Lotko, 1981a], which propagates 
upward relative to the spacecraft Furthermore, the observations allow only the two possibili- 
ties illustrated in Figure 2 The downward-pointing smaller amplitude electric field is always 
seen first, corresponding to a downward- propagating compressional mode or an upward- 
propagating rarefactive mode. In either case, the net electric field is upward, which is con- 
sistent with the particle data Lotko [1981] has shown that, in the presence of an upward ion 
beam, the downward-propagating compressional mode is damped, while the upward-propagating 
rarefactive mode can be amplified by resonant particles A rarefactive precursor to the double 
layer is also seen in the current-driven simulations discussed below. 


A train of —150 small amplitude double layers has been seen in one 45 second interval. 



One hundred fifty e4>/ T e — -1 double layers are adequate to accelerate plasma sheet electrons 
from 100 eV to 15 keV. Thus, the entire auroral potential drop might be distributed in small 
amplitude double layers 


I. Simulations of Current Driven Double Layers 

The preceding S3-3 data looks very much like the small amplitude double layers seen in 
current-driven simulations in long systems [Sato and Okuda, 1980, 1981, Hudson and Potter, 
1981; Kindel et al , 1981], There have been two distinctly different types of double layer simu- 
lations, those in which a potential is applied across the system, and those in which a current is 
imposed In the former case, whatever potential is imposed appears as a sheath at the boun- 
dary, which then detaches and propagates into the system Thus large amplitude ei f>/ T e » 1 
double layers are possible when a potential difference is imposed across a bounded system. 
When a current is imposed, current-driven instabilities result which trap ions Holes in ion 
phase space with associated negative potential spikes then appear (cf , Hudson and Potter, 1981). 
As noted by Hasegawa and Sato [1980], in the presence of a current, such a negative potential 
spike will reflect some of the electrons carrying the current, resulting in an excess of electrons 
on one side, a deficit on the other, and a net double layer potential jump. Thus it appears that 
the ion holes play a crucial role in the formation of double layers in current-driven systems 

We have examined the T e l T, dependence of double layer formation in a series of 
current-driven simulations using a system 2048 \ D long with 32000 particles for an electron 
drift equal to its thermal speed We find that holes in ion phase space and double layers form 
for TJ T, = 20 and 4 but not for T e / T, = 3 and 1 This result is related to the existence cri- 
teria for ion holes discussed below. 

II. Analytic Ion Holes and Double Layers 

We have examined time-stationary solutions to the Vlasov-Poisson equation for ion holes 
and double layers using the following distribution functions for free and trapped electrons and 
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where C= 1 A/27T and .4 = exp(— vq/2) [F(vij /2,\p / 9) + r(/3,0/0)]. /I is solved for by 
evaluating F and T, which are defined in the Appendix. The parameters «o, vo, a and (3 are 
defined below, 9 = T e / T, is the free electron-ion temperature ratio, and the potential varies 
with increasing x from 0 = 0 to 0 = -0 normalized to the free ion temperature. Schamel 
[1972] first used these distribution functions to study ion acoustic solitons. Subsequently, 
Bujarbarua and Schamel [ 1981] studied ion holes using Boltzmann electrons, which do not carry 
the current required to produce double layers in the simulations. Recently, Schamel and Bujar- 
barua [1981] have computed double layer solutions using the above distribution functions with 
current, but for a > 0, which corresponds to an ion beam vs an ion hole in the trapping 
region This type of solution permits the large amplitude e0/ T e » 1 double layers obtained 
when a potential is applied across the system in simulations and laboratory experiments. We 
will concentrate on the small amplitude (e0/ T e ~1) double layers which are connected in 
parameter space to ion hole solutions to the Vlasov-Poisson equation. These are obtained using 
the full set of distribution functions above for a < 0, which corresponds to a depletion of ions 
in the trapping region. Our primary contribution is to vary the electron drift vo, as well as 
T e / T,. Other parameters of the solutions are ion drift «o and the trapping parameters 
a — TJ T„ and f3 = TJ T et . While a and (i are formally the ratios of free to trapped particle 
temperatures, they are physically proportional to the density of ions and electrons in the trap- 
ping region. Thus, <*(/?) = 1 and mq (vq) =0 reduces the trapped distributions in (1) to 



Boltzmann, while aif. 3) > 1 corresponds to an excess and a(/3 ) <0 corresponds to a deficit of 
particles in the trapping region (cf , Figure 1, Schamel, 1972). 

Integrating Poisson’s equation once yields a simple harmonic oscillator equation for $ (x) 
in a classical potential V (0) which determines the spatial dependence of $ 

16 ( 0'(*)) 2 + Vty ) =0 ( 2 ) 

The Sagdeev potential V for ion holes and double layers is shown in Figure 3. The latter must 

satisfy V'(— 0 mm ) = 0 while the former need not This is because V'(<f>) is proportional to 

charge density, since (2) is an integral of Poisson’s equation, and there is a net charge density 

at — 0mm for a soliton On the other hand, the electric field and therefore charge density must 

vanish at — 0 mm for a laminar double layer This imposes the additional constraint 

V (-0 mm ) = 0 on double layers, which therefore occupy a subset of ion hole parameter space 

In our first study we reproduced Bujarbarua and Schamel’ s [1981] result for Boltzmann 
electrons Figure 4 is an existence diagram for ion holes in parameter space which agrees with 
their Figure 7. We present this result because our integration scheme described in the Appen- 
dix differs slightly from theirs. The dashed line indicates a cutoff where ion holes become dou- 
ble layer solutions satisfying V'(— 0 min ) = 0. This and the previous figure suggest that an ion 
hole growing in amplitude by resonant particle interaction may evolve into a double layer. 
Increasing -0 min in Figure 3a can convert the ion hole Sagdeev potential into a double layer 
Sagdeev potential like Figure 3b However, we emphasize that the calculation of Figure 4 is 
time-stationary, and the foregoing is merely a plausibility argument to suggest that when the 
ion hole amplitude becomes sufficiently large, it becomes a double layer. 

We have explored the dependence of ion hole and double layer solutions on the parame- 
ters T e / T n 1 3 = TjTg, and v 0 For /3 = 1 and v 0 = 0 corresponding to Boltzmann electrons, 
increasing T e / T, to 20 in Figure 5a and decreasing m T e f T, to 4 in Figure 5b does not change 
the available parameter space for ion hole solutions if one normalizes to e<bl T e That is, the 
solutions are always restricted to e<£/ T e < 1 for Boltzmann electrons. However, decreasing 
T e ! T, to 3 in Figure 6 changes the topology to multivalued solutions 
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We also varied the electron trapping parameter 0. In Figure 7 for /3 = 3, the parameter 
space shrinks and e(j>/ T e ~0 4 is the maximum allowed amplitude for an ion hole or a double 
layer. In Figure 8 for /3 = 0 1, the parameter space expands to allow e<f>/ T e = 2 5 For the 
Boltzmann electron case (/3 = 1), the amplitude was always restricted to e$/ T e < 1. Increas- 
ing p increases the density of electrons reflected from the flanks of the ion hole, which may 
short out the potential and reduce the maximum allowed e<j>/ T e , and conversely 

Next we examined the vo dependence. For /3 = 1 and Tj T, = 20 and 4, Figure 9a and 
Figure 9b show that the maximum allowed e0/ T e is larger for vo = 0.6 than for v 0 = 0 (see 
Figure 5a, b) The maximum allowed ion hole amplitude e4>/ T e increases with vo, but still 
remains close to unity, even for the large values of vo — 0 6 and 1 used in the current-driven 
simulations [Sato and Okuda, 1980, Hudson and Potter, 1981] As vo approaches unity 
(vo > 0 8), double layer solutions are no longer obtained at the maximum amplitude ion hole 
solution This result is discussed in the context of current-driven simulations in the following 
section For (3 = 1 and T e / T, = 3, the same kind of multivalued solutions are obtained for 
arbitrary vo as in Figure 6 

Finally, we allowed both a and /3 to assume negative values No upper limit on ion hole 
amplitude and no double layer solutions were obtained for this case 

III. Discussion 

To study the temporal evolution of an ion hole into a double layer requires examination 
of the time dependent Vlasov-Poisson equation, obtaining an evolution equation for the particle 
distribution functions and potential a' la the KdV equation Lotko [1981a] has worked on the 
nonlinear growth and damping of ion acoustic sohtons due to resonant particle interaction, and 
has shown that compressional ion acoustic sohtons can be amplified by interacting with an ion 
stream, while rarefactive ion acoustic sohtons amplify by interacting with an electron stream, 
which exists in the ion beam frame. Lotko [1981b] has shown that the downward-propagating 
compressive mode in Figure 2a will be damped because it is propagating against the ion beam. 
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while the upward propagating rarefactive mode in Figure 2b can be amplified. Whether a mode 
is amplified or damped depends on whether there are more or fewer resonant particles traveling 
faster than the wave, as in the case of linear Landau growth or damping In the saturated state 
where there is no further wave particle momentum exchange, it is possible to construct a 
potential profile which looks like Figure 2b by matching a soliton and double layer solution to 
the Vlasov-Poisson equation at the potential minimum The foregoing analysis is for ion acous- 
tic sohtons. Analysis of the temporal evolution of ion holes is in progress. 

The time-stationary analysis in the preceding section is consistent with the simulations we 
have performed. For example, the ion holes are observed to evolve into double layers with 
e<£/ T e ~\ for T e / T j = 20 and 4, but neither ion holes nor double layer formation are 
observed for TJ T, = 3 and 1, for vo = 1 initially. Although vo = 1 initially in our simula- 
tions, it drops to v 0 = 0 8 by r = 500 uj} for T e / Ti = 20, well after ion holes have formed, 
but before the double layer has formed [ Hudson and Potter, 1981] 

In summary, current-driven simulations done before solitary waves and double layers 
were seen in the S3-3 data look remarkably like the observations The simulations and time- 
stationary solutions to the Vlasov-Poisson equation suggest that these small amplitude double 
layers evolve from holes in ion phase space when the amplitude reaches e<p/ T e — 1 
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IV. Appendix 

Our treatment differs from that of Schamel and Bujarbarua [1981] in how the velocity dis- 
tributions are integrated to obtain densities, and in how those densities are integrated to obtain 
a pseudo-potential In both treatments the relevant functions are integrals over velocity space 
of the particle distribution functions multiplied by a factor exp(wo/2), where wo is the drift 
velocity of the particles Our integration of the trapped distribution is the same as that of 
Schamel [1972] and gives a function T related to the plasma dispersion function and the Daw- 
son integral 

T(x,y) = 2C exp(xy) J exp(— xt 2 /2)dt 
o 

where C= 1/V2 tt Our treatment differs in the integration of the free population Schamel 
integrates exp(v 0 2 / 2) times the free population 

/ f (v) = C exp(— '/ 2 [±(v 2 - 2<b)' h + v 0 ) 2 ] I v | > ( — 24>)' h (Al) 

over velocity space to obtain 

F(vo/2,<f>) =exp(vd/2) (1 — er/V</>) exp(c£) + K( Vq / 2, </>) (A2) 

[Schamel, 1972] where 
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K(x,y) = l/'Tir J d0Vx cos 0 exp(— ytan 2 0 + xcos 2 0) erf(Vx cos0) 

0 

We, however, make a change of variables in the velocity integral u = ± (v 2 — 20)'^ and use 

oo 

F(v 0 2 /2, 0) = 2 C f du j ~~ — rrr exp(— m 2 /2) cosh (m Vo) (A3) 

*o lu +20J" 

This can be written in two ways 

• oo 

F( v 0 2 / 2, 0) = f </« - T , - [exp (- H 2 / 2 + u Vfl) + exp (- w 2 / 2 - m vq)] (A4a) 

“o \.u l + 20 ] n 


F(vq / 2, 0) = -2 70/ 7r + Cexp(v 0 2 /2 ){ f dw{[{w + vo) 2 + 20]' /i + [(w — Vo) 2 + 20] 

o 

v 0 

x wexp(— w 2 / 2)} — 2 J* rfiv[(w— vo) 2 + 20] w wexp(— w 2 /2) (A4b) 

o 

In our evaluation, we make the change of variables « = V 2 tan0 to change the infinite 
limits to a finite interval of zero to 27r, and we numerically integrate this result. (A4a) is used 
for v 0 < 4 5, and (A4b) is used for v 0 > 4.5 

The difference in evaluation of F also leads to different forms for the integral of F used 
in the pseudo-potential. Schamel integrates (A2) over 0 for the pseudo-potential, whereas we 
integrate (A3) over 0 to find 

oo 

<7(vo /2,0) = 2cf du u [« 2 + 20] w exp(— u 2 /2) cosh (u vo) (A5) 

o 

The Sagdeev potential is then 

~0 [flexp(— V q /2) [<7( vo / 2, (0 + 0) / 6>) — G(v 0 2 /2, 0/0) + ZX0, (0 + 0)/0) 

- Z)(/3, 0/0)] + >lexp(-i/o 2 /2)[G(-Wo 2 /2,-0) - G(-u 0 2 / 2,0) + D(a,~0)]] 
where D(x,y ) == (/’(x,^) - 2-Jy/n)/ x. 

For v 0 < 4 5, (A5) is integrated just as is F in (A4a) For vo > 4 5, we use a method 
similar to that used in (A4b) 

The small amplitude expression for our F agrees with that of Schamel [1972] and the 
numerical evaluation of F agrees with the small amplitude expansion to better than 1 part in 
10 5 The numerical evaluation of G agrees with the small amplitude expansion to better than 1 
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part in 10 8 . 





Figure Captions 


Figure 1 The two perpendicular and one parallel electric field components shown Examples of double 
layers (DL), solitary waves (SW) and electrostatic ion cyclotron (EIC) waves are marked 

Figure 2 (a) Downward propagating compressive and (b) upward propagating rarefactive modes consistent 

with S3-3 electric field data / is ionosphere and M is magnetosphere 

Figure 3 Sagdeev potential V defined by (2) for (a) ion holes and (b) double layers 

Figure 4 Existence diagram for ion holes assuming Boltzmann electrons a is the ratio of free to trapped 
ion temperature and i/q is the free ion drift in the wave frame normalized to its thermal speed Negative a 
corresponds to a deficit of trapped ions (hole), see (1) The dashed line corresponds to double layer solu- 
tions 

Figure 5 (a) Same as Figure 4 for T e / T, = 20 (b) Same as Figure 4 for T e / T, — 4 

Figure 6 Same as Figure 4 for TJ T, = 3, except that solutions are now multi-valued where curves cross, 

while end points satisfy double layer criterion V (-d> m , n ) — 0 

Figure 7 Same as preceding figures with /3 > 1, which corresponds to an excess of electrons over the 
Boltzmann level in the trapping region /3 is the ratio of free to trapped electron temperature, proportional to 
the density of trapped electrons (see (1)) 

Figure 8 Same as Figure 7 with /3 < 1 , which corresponds to a deficit of electrons relative to the Boltzman 
level in the trapping region 

Figure 9 (a) Same as Figure 5a with v 0 = 0 6, where v 0 is the free electron drift in the wave frame normal- 
ized to its thermal speed (b) Same as figure 5b with vq = 0 6 
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